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Scaling of the higher-order flow harmonics: implications for 
initial-eccentricity models and the "viscous horizon" 
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(Dated: January 21, 2013) 

The scaling properties of the flow harmonics for charged hadrons v n and their ratios 
[u„/(u2) n//2 ]n>3, are studied for a broad range of transverse momenta (pr) and centrality selec- 
tions in Au+Au and Pb+Pb collisions at ^/snn = 0.2 and 2.76 TeV respectively. At relatively low 
Pt, these scaling properties are found to be compatible with the expected growth of viscous damp- 
ing for sound propagation in the plasma produced in these collisions. They also provide important 
constraints for distinguishing between the two leading models of collision eccentricities, as well as 
a route to constrain the relaxation time and make estimates for the ratio of viscosity to entropy 
density rj/s, and the "viscous horizon" or length-scale which characterizes the highest harmonic 
which survives viscous damping. 

PACS numbers: 25.75.-q, 25.75.Dw, 25.75.Ld 
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Full characterization of the transport properties of the 
strongly interacting matter produced in heavy ion colli- 
sions, is a central goal of the experimental heavy ion pro- 
grams at both the Relativistic Heavy Ion Collider (RHIC) 
and the Large Hadron Collider (LHC). Collective flow, as 
manifested by the anisotropic emission ofparticles in the 
plane transverse to the beam direction |l|, continues to 
play an important role in these ongoing efforts 
This anisotropy can be characterized, as a function of 
particle transverse momentum px and collision central- 
ity (cent) or the number of participant nucleons N part , 
by the Fourier coefficients v n ; 



dN 



n=x / 



(1) 



23 and by the pair-wise distribution in the azimuthal angle 

24 difference (A</> = <p a — </>b) between particle pairs with 

25 transverse momenta p T and p\ (respectively) [jj, LL6| ; 



dN pa 



dAcp 



x 



^2<^cos(nA0) 



(2) 



26 where 4> is the azimuthal angle of an emitted particle, and 

27 ^f n are the azimuths of the estimated participant event 

28 planes [H, [H| ; 



(cosn(<f> - * n )) 
- * m )) , n ^ m, 



(3) 



where the brackets denote averaging over particles and 
events. Here, the the starred notation is used to distin- 
guish the n-th order moments obtained relative to the ra- 
th order event plane \& m (eg. v% relative to ^2)- For flow 



driven anisotropy devoid of non flow effects, the Fourier 
coefficients obtained with Eqs. Q] and [2] are equivalent. 

Flow coefficients stem from an eccentricity-driven hy- 
drodynamic expansion of the matter in the collision zone 
a finite eccentricity e n drives uneven 
pressure gradients in- and out of the event plane, and the 
resulting expansion leads to the anisotropic emission of 
particles about this plane. The coefficients v n (pT, cent) 
(for odd and even n) are sensitive to both the initial ec- 
centricity and the specific shear viscosity rj/s {i.e. the 
ratio of shear viscosity 77 to entropy density s) of the ex- 
panding hot plasma 0-[8, 13, i3"23|- Here, it is notewor- 
thy that, for symmetric systems, the symmetry transfor- 
mation ^ up — > rp+tt, dictates that the odd harmonics 
are zero for smooth ideal eccentricity profiles. However, 
the "lumpy" transverse density distributions generated 
in individual collisions, can result in eccentricity profiles 
which have no particular symmetry, so the odd harmonics 
are not required to be zero from event to event. Fortu- 
itously, the pervasive assumption of a smooth eccentricity 
profile has hindered full exploitation of the odd harmon- 
ics until recently 

Because of the acoustic nature of anisotropic flow (i.e. 
it is driven by pressure gradients), a transparent way to 
evaluate the strength of dissipativc effects is to consider 
the attenuation of sound waves. In the presence of vis- 
cosity, sound intensity is exponentially damped e^~ r ^ Ts > 
relative to the sound attenuation length T s . This can 
be expressed in terms of a perturbation to the energy- 



momentum tensor T, 



1"' 



8T^(t) = exp 



2 V* 

Is T 



(4) 



where viscous damping reflects the dispersion relation 



2 



for sound propagation, and the spectrum of initial (t 
= 0) perturbations can be associated with the harmon- 
ics of the shape deformations and density fluctuations. 
Here, k is the wave number for these harmonics, and t 
and T are the expansion time and the temperature of 
the plasma respectively. For a collision zone of trans- 
verse size R, Eq. [4] indicates that viscous corrections for 
the eccentricity-driven flow harmonics with wavelengths 
2-KR/n for n > 1 {i.e. k ~ n/R), dampen exponentially 
as n 2 . The "viscous horizon" or length scale r v = 2nR/n v 
is also linked to the order of the highest harmonic n v 
which effectively survives viscous damping. That is, it 
separates the high frequency sound modes which are fully 
damped from those which are not (28l |. The sound hori- 
zon r s = drc s (r) , or the distance sound travels at 
speed c s (t) until flow freeze-out r/, sets the length scale 
for suppression of low frequency supcrhorizon modes with 
wavelengths 2tt Rf/n > 2r s , where Rf is the transverse 
size at sound freeze-out. Thus, the relative magnitudes of 
the flow harmonics v„ can provide important constraints 
for pinning down the magnitude of the transport coef- 
ficients rj/s and c s , as well as the "correct" model for 
eccentricity determinations 28- 30j|. 

Viscous damping for sound propagation in the plasma 
does not indicate an explicit px dependence for the rel- 
ative magnitudes of v n (cf. Eq. [4|. However, for a finite 
viscosity in the plasma, the resulting asymmetry in the 
energy-momentum tensor manifests as a correction to the 
local particle distribution (/) after freeze-out [3ll ] : 



f = fo + Sf( PT ), 



(5) 



93 where /o is the equilibrium distribution and 8f(px) is its 

94 first order correction. The latter acts as a viscous cor- 

95 rection and is known to reduce the magnitude of V2{px), 

96 especially for px > 0.7 GeV/c [3l|. The relative magni- 

97 tudes of v n (pT) are expected to be dominated by the dis- 

98 persion relation for sound propagation, albeit with some 

99 influence from 8f(j>x)- For relatively small values of rj/s, 

100 this influence on the p^-dependent viscous corrections 

101 would also be small. Thus, a characteristic scaling re- 

102 lationship between v n ,n>3(PT) and V2(px) might be ex- 

103 pected. 

104 In this letter, we investigate the scaling properties 

105 of v n (pr, cent) and the ratios [v n (PT)/{v2(PT)) n/2 ] n>3 
loe for charged hadrons produced in ultrarclativistic Au+Au 
107 and Pb+Pb collisions. We find scaling patterns that: (i) 
loa validate the viscous damping expected for sound propa- 
109 gation in the plasma created in these collisions, (ii) pro- 
no vide a constraint for distinguishing between the two lead- 
in ing eccentricity models, i.e. the Glaub er 1321 and the 

112 factorized Kharzeev-Levin-Nardi (KLN) [33M351 ] models, 

113 and (iii) point to an independent and robust method to 

114 estimate rj/s. 

us The double differential data, 7j*(pr,cent) and 
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FIG. 1. Vn/sn vs. n for charged hadrons (px ~ 1.4 GeV/c) 
produced in Au+Au collisions at ^/sjvjv = 0.2 TeV (a) and 
Pb+Pb collisions at A /s jyjy = 2.76 TeV. The v n data are taken 
from Refs. [3^1 and [37],[38| respectively for 20-30% centrality. 
The curves represent fits to the data (see text). 



117 measurements by the PHENIX collaboration, for Au+Au 
us collisions at ^/sWn 



0.2 TeV [18|, |36j, and measure- 

119 ments by the ATLAS collaboration for Pb+Pb collisions 

120 at ^J^v = 2.76 TeV [H, HI]. The Au+Au data set 

121 exploits the event plane analysis method (cf. Eq. [3]), 

122 while the Pb+Pb data set utilizes the two-particle A(f> 

123 correlation technique (cf. Eq. [2]), as well as the event 

124 plane method. Note as well that, due to partial error 

125 cancellation, the relative systematic errors for the ratios 

126 v n /{v2) n / 2 and u*/^)™/ 2 can be smaller than the ones 

127 reported for the v n values. 

128 To perform validation tests for viscous damping com- 

129 patible with sound propagation, the measured values of 

130 w n (ccnt), for each px selection, were first divided by 

131 e„(cent) and then plotted as a function of n. Monte 

132 Carlo (MC) simulations were used to compute e„(cent) 

133 from the two-dimensional profile of the density of sources 

134 in the transverse plane p s (r±), with weight w(rjj = r+" 

135 [2^|. Figs. Q] (a) and (b) show representative examples of 

136 v n lE n vs. n for charged hadrons (px ~ 1-4 GeV/c) in 

137 mid-central Au+Au and Pb+Pb collisions respectively. 

138 They confirm the exponential decrease of v n /e n with n 2 , 

139 expected for sound propagation. This "acoustic scal- 

140 ing" of v n is further confirmed by the dashed and dot- 

141 dashed curves which indicate exponential/Gaussian fits 

142 (Ae~P n ) to the data shown. 

143 Similar patterns were observed for a broad selection of 

144 centralities for px < 3 GeV/c. However, for the 0-5% 

145 and 5-10% most central Pb+Pb collisions, V2/S2 shows 

146 significant suppression relative to the empirical trend for 

147 v n vs. n, for other centralities shown by the curves in 



t (pT,cent), employed in our analysis are obtained from us Fig. [T] The fractional magnitude of this suppression is 
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FIG. 2. P vs. p T for 20-30% central Au+Au collisions at 
v /s^F = 0.2 TeV (a) and 20-30% central Pb+Pb collisions at 
y/SNN = 2.76 TeV (b). The dot-dashed and dashed curves 
indicate a 1/\/(p t ) dependence for /3 (see text). 
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149 essentially independent of pr even though «2 /£2 < V3 /£3 

150 for pr > 2 GeV/c. We interpret this suppression as an 

151 indication that, for the most central Pb+Pb collisions, 

152 the low frequency modes n < 3 exceed the superhori- 

153 zon limit, i.e. 2nRf/n > 2r s . That is, for these low 

154 frequency modes, the requirement for the maximum mo- 
ra mentum anisotropy to develop is not met, so only a frac- 

156 tion of the full anisotropy is observed. For mid-central 

157 collisions, these sound modes have shorter wavelengths 

158 which make them sub-horizon. Note that Rf = R + r s , 

159 so the order n of the low frequency modes which are sup- 

160 pressed, can serve to constrain the sound speed. 

161 For each centrality, exponential fits {Ae~P n ) to v n /e n 

162 vs. n were also made to investigate the pr-dependent vis- 

163 cous corrections attributable to 5f{pr). Figs. [2] (a) and 

164 (b) show the pr-depcndence of the j3 values extracted 

165 for 20-30% central Au+Au and Pb+Pb collisions respcc- 

166 tively; similar data trends were observed for other cen- 

167 tralitics. The dashed and dot-dashed curves in Fig. [2] 

168 show that /3 scales as 1/y/lpr) for both collision ener- 

169 gies, but the values for Pb+Pb collisions are about 25% 

170 larger. This scaling is a clear indication of the influence 



of the relaxation time [15J, |31|. Consequently it should 



173 serve as an important constraint for models. 

174 Figure [3] shows the ratios i>3/(i>2) 3 ' /2 and v^jiv-i) 2 plot- 

175 ted as a function of pr [(a) and (b)] and -/V par t [(c) and 

176 (d)] respectively, for Au+Au collisions. These ratios in- 

177 dicate an essentially flat dependence on pt, but show a 



characteristic increase with N, 



part ■ 



The same trends are 



exhibited by the Pb+Pb data (^/sWn = 2.76 TeV) with 
magnitudes comparable to those for Au+Au collisions 
for the same range of px and centrality selections. We 
interpret the fiat pr dependence of v n /{v2) n ^ 2 [for each 



FIG. 3. v 3 /(v 2 ) 3/2 vs. p T (a) and v 4 /(v 2 ) 2 vs. p T for 10-20% 
central Au+Au collisions. The bottom panels show i>3/(i>2) 3 ^ 2 
vs. iVpart (c) and ui/(d2) 2 vs. iVp ar t (d) for several pt cuts, as 
indicated. The ^2,3,4 values used for these ratios are reported 
in Ref. [3. 



183 cent] to be an indication that the py-dcpendent viscous 

184 corrections for v n are dominated by the pr-indepcndcnt 

185 contributions which stem from the dispersion relation for 

186 sound propagation. 

187 The trends for v%/(v2) 2 were found to be similar to 

188 those for v^l (v?) 2 , but the ratios v%/(v2) 2 vs. Af part are 

189 much less steep [18)]. The Af par t dependence of w|/(f2) 2 

190 and Vi/^2) 2 contrasts with the constant value of ~ 0.5 

191 predicted for perfect fluid hydrodynamics [1^, H(| , and 

192 points to the important role of the higher-order eccen- 

193 tricity moments and their fluctuations [3, [13, EH E2 . 

194 The apparent differences between ^/(i^) 2 and Vi/{v2) 2 

195 are also an indication of the important role of £4 as a 

196 driver for U4. That is, the expected contribution to V4 

197 from V2 [~ ( v 2) 2 ] does not dominate the V4 measure- 
rs ments. Note as well that V4 > v% is expected because 
199 the initial eccentricity fluctuations cause ^2 to fluctuate 

about * 4 . 

201 The flat px dependence for v n /(v2) n ^ 2 (c.f Fig. [3]) also 

202 suggests that the pr-dependent contributions to the vis- 

203 cous corrections for the ratios (v n /e n )/ (v 2 /e2) n ^ 2 essen- 

204 tially cancel, making them a reliable constraint for the ra- 

205 tios £ n /(e2)"^ 2 and consequently, an important route for 

206 distinguishing between different eccentricity models [2^] . 

207 The solid symbols in Fig. |4]show a representative set of 

208 the experimental v n j (v^) 11 ^ 2 ratios which take account of 

209 the relatively small effects of acoustic suppression. The 



4 




FIG. 4. Data comparisons to the calculated ratios (a) 
£3/(e 2 ) 3/2 vs. iV pal .t and (b) e 4 /(e 2 ) 2 vs. N part for MC- 
Glauber and MC-KLN initial geometries for Au+Au collisions 
(see text). 



210 open symbols show the corresponding eccentricity ratios 
2n obtained for the two eccentricity models. The e n values 

212 for these ratios were evaluated as described earlier. Fig. 

213 U indicates relatively good agreement between data and 

214 the enjie-if 1 !' 1 ratios, confirming the utility of v n /{v2) n ^ 2 

215 as a constraint for distinguishing between the eccentricity 
2M models (29j . 

2ia The observed scaling patterns summarized in Figs. Q] 

219 - [3] undoubtedly provide an important set of constraints 

220 for detailed comparisons to model calculations. In lieu 

221 of such calculations, we demonstrate their current util- 

222 ity for first rough estimates of the magnitude of rj/s 

223 and the viscous horizon. To this end, we employ fits to 

224 both the Au+Au and Pb+Pb data, with the fit function 



225 (A^J\pT)e V(pt) ; where the y/lpr) factors account for 

226 the influence of 5f(pr). These fits indicate that, relative 

227 to V2/S2, the magnitude of «„/£«,, n>3 decreases by more 

228 than a factor of 50 for n v ~ 7, i.e. v n /s n for n > 7 are 

229 essentially completely damped. This gives the estimate 

230 r v = 2irR/n v ~ 1.8 fm for the viscous horizon in central 

231 Au+Au and Pb+Pb collisions. 

232 The same fits allow robust extraction of the /3 values 

233 for Au+Au and Pb+Pb collisions. Here, it is noteworthy 

234 that the dependence of v„ /e n on n can provide a partic- 

235 ularly tight constraint, because it is the relative magni- 

236 tudes of v n /e n which now serve to constrain (3. These f3 

237 values have been used to extract a first rough estimate of 

238 Amj/s ^1.2 for central Au+Au collisions for (T) = 220 



MeV and t = 9 fm [12J. This rough estimate is in rea- 



240 sonable qualitative agreement with the values from prior 

241 extractions 

US M, EL Q EJ E3 • A similarly rough es- 

242 timate from the Pb+Pb data (^/snn = 2.76 TeV) gives 

243 a value for rj/s which is approximately 25% larger (cf. 

244 the larger value for /3) if we assume that the ratio T/t is 

245 roughly the same for Au+Au and Pb+Pb collisions [43J . 

246 That is, we assume that a possibly larger flow freeze- 



out time is compensated for, by a higher estimated mean 
temperature. More detailed model calculations are re- 
quired to address the values of (T) and t required for a 
more accurate estimate of rj/s. Nonetheless, our proce- 
dure clearly demonstrates the value of the relative mag- 
nitudes of v n for extraction of rj/s. 

In summary, we have presented a detailed study of the 
scaling properties of the flow coefficients v n and their ra- 
tios [v n /(v 2 ) n ^ 2 ] n >3, for Au+Au and Pb+Pb collisions 
at -y/sjvAf" = 0.2 and 2.76 TeV respectively. Within an 
empirically parametrized viscous hydrodynamics! frame- 
work, these properties can be understood to be a con- 
sequence of the acoustic nature of anisotropic flow, i.e, 
the observed viscous damping is characteristic of sound 
propagation in the plasma produced in these collisions. 
This interpretation not only provides a straightforward 
constraint for distinguishing between the two leading ec- 
centricity models, it provides a means to constrain the 
relaxation time and the sound speed, as well as to make 
independent estimates for the the averaged specific shear 
viscosity and the viscous horizon, via studies of the rela- 
tive magnitudes of v n . The observed scaling also has im- 
portant implications for accurate decomposition of flow 
and jet contributions to two-particle A4> correlation func- 
tions. This is because the higher-order harmonics can be 
expressed as a power of the high precision v% harmonic. 
It will be valuable to perform detailed viscous hydrody- 
namical model comparisons to v n and v n / (v^f 1 ! 2 for both 
identified and unidentified hadrons, as well as to establish 
the pt value which signals a breakdown of these scaling 
patterns. 
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